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Summary. The properties of X-ray emission from accreting black holes are re- 
viewed. The contemporary observational picture and current status of theoretical 
understanding of accretion and formation of X-ray radiation in the vicinity of the 
compact object are equally in the focus of this chapter. The emphasis is made pri- 
marily on common properties and trends rather than on peculiarities of individual 
objects and details of particular theoretical models. The chapter starts with discus- 
sion of the geometry of the accretion flow, spectral components in X-ray emission and 
black hole spectral states. The prospects and diagnostic potential of X-ray polarime- 
try are emphasized. Significant attention is paid to the discussion of variability of 
X-ray emission in general and of different spectral components - emission of the ac- 
cretion disk, Comptonized radiation and reflected component. Correlations between 
spectral and timing characteristics of X-ray emission are reviewed and discussed in 
the context of theoretical models. Finally, a comparison with accreting neutron stars 
is made. 



1 Introduction 

The gravitational energy of matter dissipated in the accretion flow around a compact 
object of stellar mass is primarily converted to photons of X-ray wavelengths. The 
lower limit on the characteristic temperature of the spectral energy distribution 
of the emerging radiation can be estimated assuming the most radiatively efficient 
configuration - optically thick accretion flow. Taking into account that the size of 
the emitting region is r ~ Wrg (rg is the gravitational radius) and assuming a black 

body emission spectrum one obtains: kThb ~ (ix/o"ss7rr^) ^'^^ « lAL^^ /M\1'^ 
keV. It is interesting (and well known) that Tbb scales as oc M~^^'^ . This is confirmed 
very well by the measurements of the disk emission temperature in stellar mass 
systems and around supermassive black holes in AGN. It is also illustrated, albeit 
less dramatically, by the comparison of the soft state spectra of black holes and 
neutron stars, as discussed later in this chapter. The upper end of the relevant 
temperature range is achieved in the limit of the optically thin emission. It is not 
unreasonable to link it to the virial temperature of particles near the black hole, 
kT^ir ~ GMm/r oc mc^ /{r /rg). Unlike the black body temperature this quantity 
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does not depend on the mass of the compact object, but does depend on the mass 
of the particle m. For electrons T^ir ~ 25(r/10r-g)^^ keV and it is correspondingly 
rrip/me — 1836 times higher for protons. Protons and ions are the main energy 
reservoir in the accretion flow, but for all plausible mechanisms of spectral formation 
it is the temperature of electrons that determines the spectral energy distribution of 
the emerging radiation. The latter depends on the poorly constrained efficiency of 
the energy exchange between electrons and protons in the plasma near the compact 
object. The values of the electron temperature typically derived from the spectral 
fits to the hard spectral component in accreting black holes, kTe ~ 50— 150 keV, are 
comfortably within the range defined by the two virial temperatures. However, the 
concrete value of kTe and its universality in a diverse sample of objects and broad 
range of luminosity levels still remains unexplained from first principles, similar to 
those used in the derivation of the temperature of the optically thick soft component. 

Broadly speaking, significant part of, if not the entire diversity of the spectral 
behavior observed in accreting black holes and neutron stars can be explained by 
the changes in the proportions in which the gravitational energy of the accreting 
matter is dissipated in the optically thick and optically thin parts of the accretion 
flow. The particular mechanism driving these changes is however unknown - despite 
significant progress in MHD simulations of the accretion disk achieved in recent years 
(see Chapter 10) there is no acceptable global model of accretion onto a compact 
object. To flnish this introductory note, I will mention that non-thermal processes 
in optically thin media (e.g. Comptonization on the non-thermal tail of the electron 
distribution) may also contribute to the X-ray emission from black holes in some 
spectral states. 
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Fig. 1. Energy spectra of a black hole Cyg X-1 (left) and a neutron star 4U 1705-44 
(right) in the soft and hard spectral states (left panel adopted from [19]). 



X-ray emission from black-hole binaries 



3 



2 Geometry and spectral components 

The contributions of optically thick and optically thin emission mechanisms can be 
easily identified in the observed spectra of accreting black holes as soft and hard 
spectral components (Fig. [T|. Depending on the spectral state of the source one of 
these components may dominate the spectrum or they can coexist giving compara- 
ble contribution to the total emission. The soft component is believed to originate 
in the geometrically thin and optically thick accretion disk of the Shakura-Sunyaev 
type . To zero-th order approximation, its spectrum can be regarded as a super- 
position of black body spectra with the temperatures and emitting areas distributed 
according to the energy release and the temperature profile of the accretion disk. 
The simplest example of this is the so-called "multicolor disk blackbody model" in- 
troduced by ^5]. Although valid for one particular inner torque boundary condition 
allowing easy integration of the total fiux, this model has been widely and efficiently 
used in the era of more limited computer resources due to its simplicity, speed and 
early integration in the XSPEC spectral fitting package. To achieve a higher degree 
of accuracy one would need to consider more realistic inner boundary conditions, 
deviations of the gravitational potential from the Newtonian and to account for such 
effects as distortion of the black body spectrum due to Thomson scatterings in the 
upper layers and in the atmosphere of the disk, Doppler effect due to rotation of 
the matter in the disk etc. A number of models has been proposed to include these 
effects [121 1501 H71 FT], many of them currently implemented in the XSPEC package. 

The inadequacy of the optically thick emission mechanism in describing the 
hard spectral component present in both spectra in Fig. [1] can be easily demon- 
strated. Indeed, the size of the emitting region required to achieve a luminosity 
of \og{Lx) ~ 37 — 38 with black body temperature of ~ several tens of keV is 
Rem ~ {Lx/nasBT^y^'^ ~ meters (T30 = Ttt/SO keV). Obviously this 

is much smaller than the size of the region of main energy release near an accret- 
ing black hole, ~ 50rg ~ 500 — 1500 km. The bremsstrahlung emission from a ~ 
uniform cloud of hot plasma with a filling factor close to unity can not deliver the 
required luminosity either, because the cx density dependence of its emissivity 
makes it very inefficient in the low density regime characteristic of the optically thin 
accretion flow. Indeed, the optical depth of such a plasma cloud of size ~ 10— lOOr^ 
and emission measure cm required to explain observed hard X-ray 

luminosity of ~ lO^'^ erg/s would greatly exceed unity. However, bremsstrahlung 
emission may play a role in the advection dominated accretion flow [38] in the low 
M regime of quiescent state of accreting black holes {log{Lx) < 32 — 33). 

It has been understood early enough that Comptonization is the most plausible 
process of formation of the hard spectral component [551 156]. Thanks to its linear 
dependence on the gas density, Comptonization of soft photons on hot electrons 
in the vicinity of the compact object can efficiently radiate away the energy dissi- 
pated in the optically thin accretion flow and successfully explain the luminosity 
and overall spectral energy distribution observed in the hard X-ray band. Moreover, 
Comptonization models of varying degree of complexity do satisfactorily describe 
the observed broad band energy spectra of black holes to the finest detail (e.g. [65]). 
It is especially remarkable given the high quality of the X-ray data which became 
available in the last decade from observations of recent and current X-ray observa- 
tories, such as Compton GRO, RossiXTE and Chandra, XMM and Swift. 
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Fig. 2. The three main components of the X-ray emission from an accreting black 
hole (top) and a plausible geometry of the accretion flow in the hard spectral state 
(bottom). 



The Comptonization site - cloud of hot (thermal or non-thermal) electrons is 
often referred to as a "corona". Although it is generally accepted that the Comp- 
tonizing corona has to be located in the close vicinity of the compact object, there 
is currently no broad consensus on the detailed geometry of the region. The numer- 
ical simulations have not reached the degree of sophistication required to perform 
full self-consistent global simulations of the accretion flow. Two physically motivated 
geometries are commonly considered: (i) the "sombrero" configuration (to my knowl- 
edge first introduced in 41 ) and (ii) the patchy, flaring corona above the accretion 
disk. 

One of the variants of the "sombrero" configuration is depicted in FigH] It is as- 
sumed that outside some truncation radius the accretion takes place predominantly 
via the optically thick and geometrically thin accretion disk, whereas closer to the 



X-ray emission from blaclt-hole binaries 



5 



compact object the accretion disk is transformed into a hot optically thin and ge- 
ometrically thick flow with the aspect ratio of H/R ~ 0.5 — 1. The soft (optically 
thick) and hard (optically thin) spectral components are formed in the accretion 
disk and the hot inner flow correspondingly. The value of the truncation radius can 
be inferred from the observations. Although their interpretation is not unique and 
unambiguous, the plausible range of values is between ~ 3 and a few hundred grav- 
itational radii. There is no commonly accepted mechanism of truncation of the disk 
and formation of the corona, with a number of plausible scenarios having been inves- 
tigated recently. Among the more promising ones is the evaporation of the accretion 
disk under the effect of the heat conduction. It was initially suggested to explain 
quiescent X-ray emission from cataclysmic variables i32j and was later applied to 
the case of accretion onto black holes and neutron stars [33]. It not only provides a 
physically motivated picture describing formation of the corona and destruction of 
the optically thick disk but also correctly predicts the ordering of spectral states vs. 
the mass accretion rate. Namely, it explains the fact that hard spectra indicating 
prevalence of the hot optically thin flow are associated with lower M values, whereas 
the optically thick disk appears to dominate the photon production in the accretion 
flow at higher M. As a historical side note I mention that in the earlier years of 
X-ray astronomy the presence of the hot optically thin plasma in the vicinity of the 
compact object was often associated with disk instabilities, therefore such a behav- 
ior appeared puzzling to many astrophysicists in the view of the M dependence of 
these instabilities. 

Another geometrical configuration considered in the context of hard X-ray emis- 
sion from black holes is a non-stationary and non-uniform (patchy) corona above 
the optically thick accretion disk. This scenario has been largely inspired by the 
suggestion by Galeev, Rosner and Vaiana [14], that a magnetic field amplified in the 
hot inner disk by turbulence and differential rotation may reach the equipartition 
value and emerge from the disk in the form of buoyant loop-like structures of solar 
type above its surface. These structures may lead to the formation of a hot mag- 
netically confined structured corona similar to the solar corona, which may produce 
hard emission via inverse Compton and bremsstrahlung mechanisms. The model 
could also explain the faintness of the hard emission in the soft state as a result of 
efficient cooling of plasma in the magnetic loops via inverse Compton effect due to 
increased flux of soft photons at higher M . 

Although the original Galeev, Rosner and Vaiana paper was focused on the 
thermal emission from hot plasma confined in buoyant magnetic loops, the latter are 
a plausible site of particle acceleration responsible for the non-thermal component in 
the electron distribution. This is expected on theoretical grounds and is illustrated 
very well observationally by the presence of the non-thermal emission component 
in the spectra of solar flares (see e.g. [TO] and references therein). On the other 
hand, a power law -shaped hard X-ray component is commonly observed in the 
spectra of black hole candidates (Fig(3]l . It reveals itself most graphically in the soft 
spectral state but may be as well present in the hard state, along with the thermal 
Comptonized spectrum (e.g. [34]). This power law component has a photon index of 
7^ ~ 2 — 3 and, although it may extend into the several hundred keV - MeV range 
(e.g. [55]), it is relatively unimportant energetically, contributing a small fraction 
to the total radiation output of the black hole. This is in contrast with the hard 
spectral component due to thermal Comptonization, observed in the hard spectral 
state which accounts for the most of the source luminosity. 
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Fig. 3. Broad band spectra of Cyg X-1 in the soft and hard spectral state based on 
the data of BeppoSAX and Compton GRO missions. Adopted from [34j . 




Fig. 4. Comptonization spectrum in the case of large temperature and small optical 
depth. Left panel: kTbb = 100 eV, kT, = 300 keV, r = 0.01,0.05,0.1 from the 
bottom to the top; right panel: kTbt = 500 eV, fcTe = 1 MeV, r = 0.001,0.01,0.1. 
For comparison, the top most spectrum on the left panel shows case of r = 1. The 
left-most peak in all spectra is made of seed photons which left the Comptonization 
region without scatterings. The dashed line is a power law with a photon index of 
r = 3. These spectra are the results of Monte-Carlo simulations assuming spherical 
geometry. 
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There have been attempts to explain the power law component as a result of 
Comptonization in a media with Maxwellian distribution of electrons of large tem- 
perature and small optical depth. This seems to be unlikely as in this case humps 
due to individual scattering will be seen in the output spectrum, in contrast with 
the smooth power law spectra observed. This qualitative consideration is illustrated 
by results of Monte-Carlo simulations shown in Fig|4l 

The sombrero configuration is often associated with a predominantly thermal 
distribution of electrons, whereas the solar-like flares above the accretion disk may 
be the site of electron acceleration producing non-thermal electron distributions and 
power law-like Comptonization spectra. Thermal Comptonization is believed to be 
the main mechanism in the hard spectral state, whereas non-thermal Comptoniza- 
tion is probaly relevant in the soft state (Fig[3]). On the other hand, observations 
often indicate presence of both thermal and non-thermal Comptonization compo- 
nents suggesting that the two types of corona may co-exist (e.g. [24]). 

Due to heating of the disk by Comptonized radiation from the corona and soft 
photon feedback, the uniform thermal corona above the optically thick accretion 
disk can not explain observed hard spectra with photon index of _r ~ 1.5 — 2.0. It 
has been first proposed by Sunyaev and Titarchuk [57] that presence of cool opti- 
cally thick media - e.g. accretion disk or surface of the neutron star, in the vicinity 
of the Comptonization region will affect the parameters of the latter and, conse- 
quently, the shape of the outgoing Comptonized radiation. Indeed, some fraction 
of the Comptonized radiation will be returned to the accretion disk increasing its 
temperature and, consequently, the soft photon flux to the Comptonization region. 
This in turn will increase the cooling rate and will decrease the electron temperature 
in the Comptonization region leading to softer and steeper spectra. I will further 
illustrate it by the following simple quantitative consideration. In the sandwich-like 
geometry, assuming moderate Thompson optical depth of the corona, tt ~ 1, a frac- 
tion / ~ 1/2 of the Comptonized emission will be returned to the accretion disk. Of 
this, a fraction of 1 — a (a ~ 0.2 - albedo) will be absorbed and will contribute to the 
heating of the accretion disk, adding to its heating due to the gravitational energy 
release. Ignoring the latter, the luminosity enhancement factor in the Comptoniza- 
tion region, defined as a ratio of its total luminosity to the luminosity of the seed 
photons (see more detailed discussion in section ing will he A~ {l — a)~^f~^ < 2.5. 
As well known [IT], the luminosity enhancement factor is intimately related to the 
Comptonization parameter y and the photon index F of the Comptonized radiation. 
The above constrain on A implies 7^ > 2.3 (see Fig|15p. which is steeper than the 
hard state spectra typically observed in black holes. This conclusion is confirmed by 
the full treatment of the Comptonization problem in "sandwich" geometry 23'. In 
order to produce a harder Comptonized spectrum, the value of the feedback coefB- 
cient / needs to be reduced. This is achieved, for example, in a non-uniform, patchy 
and/or non-stationary corona. Another example suggested by ^6 involves bulk mo- 
tion of the corona with mildly relativistic velocity away from the disk reducing the 
feedback coefficient due to the aberration effect (section |8.2|I . A uniform stationary 
corona above the accretion disk can, in principle, be responsible for the steep power 
law component often detected in the soft state, although a non-thermal electron 
distribution may be a more plausible explanation, as discussed above. To conclude, 
I also note that considerations of a similar kind involving the neutron star surface 
can explain the fact that the neutron star spectra are typically softer than those of 
black holes (e.g. [57). section O. 
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3 Spectral states and geometry 

The existence of different spectral states is a distinct feature of accreting X-ray 
sources, independently of the nature of the compact object fFig |ll3|) . Affhough their 
phenomenology is far richer, for the purpose of this chapter I will restrict myself to 
the simple dichotomy between soft (high) and hard (low) spectral states and refer 
to the next chapter of this book for a more detailed discussion. As no global self- 
consistent theory/model of accretion exists, all theories explaining spectral states 
have to retreat to qualitative considerations. These considerations, although phe- 
nomenological in nature, usually are based on numerous observations of black hole 
systems, simple theoretical arguments and some simplified solutions and simulations 
of the accretion problem. Described below is a plausible, although neither unique nor 
unanimously accepted, scenario of this kind based on the "sombrero" geometry of 
the accretion flow. There is a number of cartoons and geometry sketches, illustrating 
this and other scenarios which I will not repeat here and will refer the interested 
reader to original work, e.g. [69] 

As obvious from Fig[Tl the spectral states phenomenon is related to the redistri- 
bution of the energy released in the optically thick and optically thin components of 
the accretion flow. In the sombrero configuration one may associate spectral state 
transitions with change of the disk truncation radius - the boundary between the 
outer optically thick accretion disk and the inner optically thin hot flow. 

In the soft (aka disk-dominated) spectral state, the optically thick accretion 
disk extends close to the compact object, possibly to the last marginally stable 
Keplerian orbit (r = Svg for a Schwarzschild black hole), leaving no "room" for 
the hot optically thin flow. Therefore the major fraction of the accretion energy is 
emitted in the optically thick accretion disk giving rise to a soft spectrum of the 
multicolor blackbody type. The magnetic activity at the disk surface may (or may 
not) produce a hard power-law like tail due to non-thermal Comptonization in the 
corona. As discussed in the previous section, this power law component has a steep 
slope _r ~ 2 — 3 and is relatively insigniflcant energetically. 

In the hard (aka corona-dominated) spectral state, the accretion disk truncates 
at a distance of ~ 50— lOOrg or further from the compact object. The major fraction 
of the gravitational energy is released in the hot inner flow. Comptonization of the 
soft photons emitted by the accretion disk on the hot thermal electrons of the in- 
ner flow leads to the formation of the hard spectrum of the shape characteristic for 
unsaturated thermal Comptonization. The typical parameters in the Comptoniza- 
tion region - hot inner flow, inferred from observations are: electron temperature of 
Te ~ 100 keV and Thompson optical depth of Te ~ 1. The significance of the soft 
backbody-like emission from the optically thick disk as well as of the non-thermal 
emission due to magnetic flares at its surface varies depending on the disk truncation 
radius, increasing as the disk moves inwards. There is evidence that both thermal 
and non-thermal hard components may co-exist in the hard state in the certain 
range of the disk truncation radii, as suggested in [24] . 

4 Reflected emission 

After escaping the corona, a fraction of the Comptonized photons may be inter- 
cepted by the optically thick accretion disk. Part of the intercepted radiation will 
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be dissipated in the disk, via inverse Compton-efltect and photoabsorption on heavy 
elements, contributing to its energy balance, the remaining part will be reflected due 
to Compton scatterings 0]. The disk albedo depends on the photon wavelength and 
on the chemical composition and ionization state of the disk material. Expressed in 
terms of energy flux it is ~ 0.1 — 0.2 in the case of neutral matter of solar abundance 
[31) . Some fraction of the photobsorbed emission will be re-emitted in the fluores- 
cent lines of heavy elements and may also escape the disk. Thus, combined effects of 
photoabsorption, fluorescence and Compton scattering form the complex spectrum 
of the reflected emission, consisting of a number of fluorescent lines and K-edges of 
cosmically abundant elements superimposed on the broad Compton reflection hump 
[H [15] (Fig. [51 . The peak energy of the latter depends on the shape of the incident 
continuum and for a typical spectrum of a black hole in the hard state is located at 
- 30 keV (Fig.[2|. 




Energy, keV 

Fig. 5. Spectrum of emission reflected from an optically thick slab of neutral matter 
of solar abundance. The solid line shows the spectrum of the incident power law 
emission. 

Due to the dependence of the photoabsorption cross-section, element abundances 
and fluorescent yield on the atomic number, the strongest among the spectral fea- 
tures associated with heavy elements are the K-edge and fluorescent K-a line of 
iron. For reflection from the neutral media, the centroid of the iron line is at 6.4 keV 
and its equivalent width computed with respect to the pure reflected continuum is 
~ 1 keV. For a distant observer, the reflected emission is diluted with the primary 
Comptonized radiation and thermal emission of the accretion disk, adding complex- 
ity to the observed spectra. It leads to the appearance of characteristic reflection 
features in the spectra of X-ray binaries (Fig. [5} - the fluorescent line from iron, 
iron K-edge and a broad Compton reflection hump at higher energies [4] [15] . Their 
shape is further modified, depending on the ionization state of the matter in the disk 
[48) . by strong gravity effects and Doppler and aberration effects due to the Keple- 
rian rotation of the disk e.g. [13) . The observed line usually has non-zero intrinsic 
width and the K-edge is never sharp (smeared edge). Their amplitude depends on 
the ionization state and on the solid angle of the reflector as seen from the source 
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of the primary radiation. This can change with the spectral state of the source. The 
typical value of the iron line equivalent width is ~ 50 — 300 eV. 

The presence of theses features makes the spectra deviate from a power law 
shape, expected for Comptonized radiation with parameters typical for black holes. 
As their amplitude depends on the relative configurations of the corona and the 
accretion disk and physical conditions in the disk, they have great diagnostics po- 
tential for studying the accretion geometry in different spectral states. This is further 
discussed in Sect. |8] 



5 Polarization of X-ray emission 

As the distribution of matter in the accretion flow is not spherically symmetric, one 
may expect some degree of polarization of emerging X-ray radiation. In the conven- 
tional scenario of formation of X-ray radiation in the vicinity of a compact object the 
polarization is caused by Thompson scatterings of photons on free electrons in the 
accretion disk and hot corona and is predicted to be present at the moderate level 
of ~ several percent in all three main components of the X-ray emission: thermal 
emission of the optically thick disk, Comptonized emission and reflected component. 
A considerably larger degree of polarization may be expected in some alternative 
scenarios, for example in some versions of the jet scenario, discussed later in this 
book. 

The degree of polarization of the Comptonized emission depends strongly on the 
geometry of the corona, the location of the sources of the soft seed radiation and 
the viewing angle [581 164) . As calculated by Sunyaev & Titarchuk 58 , for the case 
of Comptonization in the disk it lies in the range between and ~ 12%, depending 
on the optical depth of the corona and viewing angle, i.e. it can slightly exceed the 
maximum value of 11.7% for a pure-scattering semi-infinite atmosphere. Although 
the disk geometry of the corona is unlikely, this result suggests that a moderate 
degree of polarization of the Comptonized continuum, of the order of ~few per cent, 
may be expected. 

Obviously, polarization should also be expected for the reflected emission from 
the accretion disk. The pure reflected component can be polarized to ~ 30% [30) . 
The degree of polarization drops to ~ few per cent when it is diluted with the 
Comptonized radiation and the thermal emission from the disk. The degree of po- 
larization is a strong function of photon energy - it is low at low energy < 10 keV 
and reaches its maximum of ~ 5% (for a 60 degree viewing angle) at ~ 30 — 50 keV, 
where the contribution of the Compton-reflected continuum to the overall spectrum 
is maximal [30) . 

The thermal emission generated inside the optically thick accretion disk is ini- 
tially unpolarized but attains polarization as a result of Thompson scatterings in 
the atmosphere of the accretion disk 9, 28 . Calculations 28 show that the degree 
of polarization varies from 0% (face-on disk) to ~ 5% (edge-on). The polarization 
degree depends on the photon energy, being largest in the Wien part of the thermal 
spectrum, where scatterings play the most important role. A remarkable property of 
polarized radiation from the vicinity of a compact object is that due to relativistic 
effects the polarization angle also becomes dependent on the photon energy [91 128). 
This is in contrast with the classical approximation, where symmetry considerations 
require that the polarization direction is coaligned with the minor or major axis 
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of the disk projection on the plane of the sky. The amphtude and shape of these 
dependences are sensitive to the disk inchnation and the spin of the black hole. In 
combination with spectral information, this can be used to resolve the degeneracy 
between the black mass, spin and disk inclination [28] . 

Polarization measurements are yet an unexplored area of high energy astro- 
physics and have great diagnostics potential. With the new generation of polarimet- 
ric detectors, X-ray polarimetry will become a powerful tool to study the geometry 
of the accretion flow and the properties of the compact object in accreting systems. 

6 Variability 

Variability of X-ray emission is a common and well-known property of X-ray bina- 
ries. The amplitude of flux variations depends on the time scale and photon energy 
and can be as large as ~ 20 — 30% fractional rms. The dependence of the variability 
amplitude on time scale is conventionally characterized by the power density spec- 
trum, which is a square of the Fourier-amplitudes of the time series, renormalized in 
order to give the answer in desired units, for example (fractional rms)^/Hz. Power 
density spectra of X-ray sources often reveal a number of rather narrow features of 
various widths superimposed on a broad band continuum of aperiodic variations, 
suggesting that both resonances of various degree of coherence as well as stochastic 
processes of aperiodic nature contribute to the observed flux variations. Some of the 
narrow features can be clearly associated with the spin frequency of the neutron 
star or orbital or precession frequency of the binary. The nature of others, broadly 
reffered to as quasi-periodic oscillations (QPO), is poorly understood. The models 
range from the beat frequency model proposed by A.Alpar and J.Shaham soon 
after the discovery of the first QPOs in neutron stars, and the relativistic preces- 
sion model [54] relating QPO frequencies with fundamental frequencies of Keplerian 
orbits in strong gravity to models employing global oscillation modes in the accre- 
tion disk [l] [61] . QPOs were first discovered in neutron star systems, but are also 
commonly observed in accreting black holes. Their rich phenomenology is well doc- 
umented elsewhere (e.g. 62 ) and is not discussed in this chapter, which will focus 
on the aperiodic variability continuum. 

I will conclude these introductory remarks with the following side note. The 
interpretation of the energy spectra of celestial X-ray sources has been greatly facil- 
itated by the fact that a number of simple physical concepts could be employed in 
a straightforward manner, such as black-body emission, Comptonization, photoab- 
sorption and fluorescence etc. No equivalent concepts are easily available to help in 
interpreting the power density spectra. This may explain why our understanding of 
the power density spectra lags significantly behind our understanding of the energy 
spectra. Indeed, spectral analysis can rely on a number of physically motivated and 
elaborate models which successfully describe high quality data from modern obser- 
vatories. The interpretation of the power density spectra, on the contrary, has just 
started to advance beyond "numerology" and simple association of QPO frequen- 
cies with characteristic frequencies of test particles in the (strong) gravitational field. 
However, it is obvious that timing information (along with the polarization measure- 
ments) presents a completely different dimension which has to be taken into account 
in validating any model of formation of radiation in the vicinity of a compact object. 
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Fig. 6. Broad-band power density spectra of Cyg X-1 and Cyg X-2. Tiie spectra 
are obtained from non-simultaneous RXTE/ASM, EXOSAT/ME and RXTE/PCA 
data. In the left panel, the open circles at log(/) > —4 (red in the color version of 
the plot) show the 1996 high state PDS, the filled circles (blue) are hard spectral 
state PDS. The solid lines extending through entire plots in the left and right panels 
show a. Pv (X u'^'^ power law; normalization is the same on both panels. The second 
power law in the left panel giving a better approximation to the soft state data of 
Cyg X-1 is oc u'^ power law. The very high frequency data log(/) > 1 are from 



6.1 Propagation of M fluctuations in the accretion disk 

The remarkable characteristics of aperiodic variability is its breadth in the frequency 
domain. As illustrated by the power density spectra of Cyg X-1 and Cyg X-2 shown 
in Fig. |S] flux variations in the frequency range from log / < —8 to log / ~ 2 — 3 are 
present. This suggests that variations of the mass accretion rate M on an extremely 
broad range of time scales are present in the innermost region of the accretion 
flow, r < 50 — 100 rg, where X-ray emission is produced. It is to be compared 
with characteristic time scales in this region which are limited by two extremes 
- the dynamical time scale that is of the order of the Keplerian orbital time tx ~ 
0.3 (M/10Mq) {r /50rg)^^^ sec and the viscous time scale tvisc ~ {h/r)~^ Q~^^ = 
(2-Ka)~^ {h/r)~^ tK {ce is viscosity parameter, h - disk thickness). The latter is in 
the range between tyisc ~ 10'* ~ 10"^ sec for a gas pressure dominated Shakura- 
Sunyaev disk {h/r ~ 10~^) and tyisc ~ 10 tj? ~ 1 — 10 sec for a thicker hot flow with 
the aspect ratio of h/r ~ 1/3. 

From the point of view of characteristic time scales, the high-frequency varia- 
tions could potentially be produced in the vicinity of the compact object. Longer 
time scales, on the contrary, exceed by many orders of magnitude the longest time 
scales in the region of the main energy release and can not be generated there. The 
low frequency M variations have to be generated in the outer parts of the accretion 
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flow and be propagated to the region of the main energy release where they are con- 
verted into variations of the X-ray flux. The power spectra shown in Fig. |6] maintain 
the same power law shape over a broad frequency range, suggesting that the same 
physical mechanism is responsible for flux variations at all frequencies. The same 
slope and normalization of the power law component in the power density spectra of 
different sources - black holes and neutron stars - suggests that this mechanism is a 
property of the accretion disk and does not depend on the nature of the compact ob- 
ject. A plausible candidate for such mechanism may be viscosity fluctuations caused 
for example by MHD turbulence in the accretion disk, as suggested in [27]. Viscosity 
fluctuations lead to variations of the mass accretion rate which after propagating 
into the innermost region of the accretion flow are transformed into variations of 
the X-ray flux. 

However, because of the diffusive nature of disc accretion 05], fluctuations 
of mass accretion rate on time scales much shorter than the diffusion time scale 
will be not propagated inwards, but will be damped close to the radius at which 
they originated. As demonstrated in [8], the amplitude of the fluctuations on 
the time scaler will be significantly suppressed at the characteristic length scale 
At jr ~ ^ r/tvisc- As the viscous time scale depends quadratically on the disk thick- 
ness, tvisc ~ {h/r)~^ ^K^j fluctuations on the dynamical t ^ td ^ or thermal 
T ^ tth ^ (ai?^)"^ time scales will be damped in the thin disk with h/r <C 1 after 
traveling a small distance in the radial direction, Ar/r ^ h/r, and will never reach 
the region of the main energy release. The inner region itself can not generate M 
fluctuations on the dynamical or thermal time scales either. Indeed, the coherence 
length for fluctuations on a time scale r ~ Qj} is Ar/r ~ h/r ~ 10~^ for the gas 
pressure supported Shakura-Sunyaev disk. Therefore N ~ r /Ar ~ 100 independent 
annuli regions will contribute to the observed flux, their variations being uncorre- 
lated. The contribution of each region to the total flux is small, ~ lO"'^, and alone 
can not cause significant variability of the total fiux. Furthermore, in their combined 
emission uncorrelated fiuctuations will be suppressed by a factor of ~ l/\/7V ~ 10~^ 
due to the averaging effect. Thus, for the geometrically thin disc, fiuctuations of vis- 
cosity (or mass accretion rate) on the dynamical or thermal time-scales will not 
contribute significantly to observed variability of the X-ray flux. In order to lead to 
significant modulation of the X-ray flux, viscosity or accretion rate fluctuations have 
to be propagated inwards significant radial distance and cause fiuctuations of mass 
accretion rate in a significant range of smaller radii, including the region of the main 
energy release. This is only possible for fiuctuations on time scales equal or longer 
than the viscous time of the disk at the radius, where fiuctuations are "inserted" 
into the accretion flow [271 [8]. 

Thus, the standard Shakura-Sunyaev disk plays the role of a low-pass filter, at 
any given radius r suppressing M variations on time scales shorter than the local 
viscous time tvisc('"). The viscous time is, on the other hand, the longest time scale 
of the accretion flow, therefore no significant fiuctuations at longer time scales can 
be produced at any given radius. Hence, a radius r contributes to M and X-ray fiux 
variations predominantly at a frequency / ~ t^\ac{r)~^ [271 [8]. The broad range of 
variability time scales observed in the X-ray power density spectra (Fig[6]) is ex- 
plained by the broad range of radii at which viscosity fiuctuations are generated. 
If viscosity fiuctuations at all radii have the same relative amplitude, a power law 
spectrum P^, oc will naturally appear [27], in qualitative agreement with obser- 
vations (Fig[S]). The picture of inward-propagating fiuctuations outlined above also 
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successfully explains the observed linear relation between rms amplitude of aperi- 
odic variability and total X-ray flux in black holes [63) and the nearly logarithmic 
dependence of the time lag between time series in different energy bands on the 
photon energy ^5] , 

6.2 Very low frequency break and accretion disk corona 

Owing to the finite size of the accretion disk, the longest time scale in the disk is 
restricted by the viscous time on its outer boundary tvisc{Rd)- Below this frequency, 
X-ray flux variations are uncorrelated, therefore the power density spectrum should 
become flat at / < tvisc{Rd)~^ ■ This explains the low frequency break clearly seen 
in the power spectrum of Cyg X-2 (Fig. [6} . If there are several components in the 
accretion flow, for example, a geometrically thin disk and a diffuse corona above it, 
several breaks can appear in the power spectrum at frequencies corresponding to 
the inverse viscous time scale of each component. It can be easily demonstrated [22] 
that the break frequency is related to the orbital frequency of the binary via: 

fbreak = 37ra(l + q)~^^^ (h/rf {Ra/a)~^^^ fort 

For low-mass X-ray binaries (Roche-lobe fliling systems) this becomes fbreak oc 
(h/r)^ forb- Such very low frequency breaks indeed are observed in a statistically 
representative sample of low mass X-ray binaries, and the break frequency is in- 
versely proportional to the orbital period of the binary (Fig. [3 [22]). 
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Fig. 7. Relation between the break frequency (~ t~l^^ at the outer edge of the accre- 
tion disk) and the orbital frequency of the binary system. The shaded area shows the 
dependence /vise vs. /orb expected for the standard gas pressure supported Shakura- 
Sunyaev disk (h/r ~ few x 10~^). Straight solid and dashed lines are predictions 
for larger values of the disk thickness h/r, as indicated by the numbers on the plot. 
The two points for Cir X-1 correspond to the original data (open circle) and after 
correction for the eccentricity of the binary orbit (fiUed circle). Adopted from [22] . 
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However, measured values of the break frequency imply that the viscous time of 
the accretion flow is a factor of > 10 shorter than predicted by the standard theory 
of accretion disks (cf. shaded area in Fig. [Tjl . This suggests that significant fraction 
of the accretion M occurs through the geometrically thicker coronal flow above the 
standard thin disk. Note that the existence of the Shakura-Sunyaev disk underneath 
the coronal flow is required by optical and UV observations of LMXBs indicating 
the presence of the optically thick media of extent comparable with the Roche-lobe 
size [22] . The aspect ratio of the coronal flow implied by the fbr/ forb measurements, 
h/r ~ 0.1, corresponds to a gas temperature of T ~ 0.01 Tvir- The corona has 
moderate optical depth in the radial direction, tt ~ 1, and contains < 10% of the 
total mass of the accreting matter (but the fraction of M is much larger, probably 
~ 0.5). These estimates of temperature and density of the corona are in quantitative 
agreement with the parameters inferred by the X-ray spectroscopic observations by 
Chandra and XMM-Newton of complex absorption/emission features in LMXBs 
with large inclination angle. 

In this picture, the red noise (power law) component of the observed variability 
of the X-ray flux is defined by the viscosity and/or M fluctuations generated in the 
coronal flow rather than in the geometrically thin disk. Fluctuations produced in 
the standard Shakura-Sunyaev disk have much smaller amplitude due to a ~ 10 
times longer viscous time scale (see discussion in Sect. I6.1|l and do not contribute 
significantly to the observed variability of X-ray flux. This is further supported by 
the lack of variability in the soft black body component in Cyg X-1, as discussed 
below. 

6.3 High and very high frequencies 

In addition to the power law component of red noise dominating in the larger part 
of the low frequency range, power density spectra usually have significant excess of 
power in the high frequency end, / > 10~^ — 10~^ Hz (Fig.[Sl). This excess power is 
often referred to as a "band-limited noise", due its limited scope in frequency. It is 
this part of the power spectrum that is the subject of investigation in the majority 
of "standard" timing analysis projects, for example in the ones based on a typical 
RXTE/PCA observation. 

In black holes, this component is most prominent in the hard spectral state and 
usually disappears in the soft state (Fig. [HI [HI)- Plotted in the vP,^ units, it shows 
several broad humps, separated by a ~ 0.5 — 1 decade in frequency (Fig. [8]). The 
centroid frequencies of these bumps vary with time but they usually do correlate 
with each other [171 166] . It has been proposed that they may be identified with the 
(much more narrow) QPO peaks observed in the power spectra of accreting neutron 
stars. When the proper identifications are made, black hole and neutron star systems 
appear to follow the same global relation between QPO frequencies [421 1661 15^ [5]. 
This suggests that the features on the power density spectra and their frequencies 
are the property of the accretion disk, their existence being unaffected by the nature 
of the compact object. 

Based on the amplitude and timescale considerations, it is plausible to link the 
origin of the band-limited noise component to the inner hot optically thin flow in 
the region of the main energy release. For a IOMq black hole, the characteristic time 
scales in this region {r ~ lOOrg) are: ~ 1 sec and t^isc ~ 10 sec, the latter was 
computed assuming a ~ 0.2 and h/r ~ 0.2. This suggest that the two main peaks 
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Fig. 8. Power density spectra of Cyg X-1 in tlie soft and iiard spectral states. The 
right hand panel shows the very high frequency part of the power density spectrum. 
The power is shown in units of uP^ (left) and u^Pi, (right). Adopted from [19] and 
[44] respectively. 



observed in the i/Pv spectra may (or may not) be related to the dynamical and 
viscous time scales at the inner edge of the truncated accretion disk and their origin 
may be linked to the interaction between the geometrically thin outer accretion disk 
and optically thin hot inner flow. Note that what appears to be the lower frequency 
hump in the uPi, plot is in fact a break, below which the power density distribution is 
flat, down to the frequency where the red noise component becomes dominant (Fig. 
[5|. As a flat power density distribution means lack of correlation between events on 
the corresponding time scales, we can draw an analogy with the very low frequency 
breaks associated with the viscous time scale at the outer edge of the disk (Sect. 
I6.2|l . This further supports the association of the low frequency hump in the vP^ 
plot. Fig. [8| (= break in the power density spectrum) with the viscous time scale at 
the outer boundary of the hot inner flow. One can estimate the truncation radius 
of the geometrically thin disk (= outer radius of the hot inner flow in the sombrero 
configuration. Sect. [2]) equating the viscous time with the inverse break frequency, 
/j^^ ~ 5 sec: rtr ~ lO^^g assuming a = 0.2 and the aspect ratio of the hot inner 
flow h/r — 1/3. This number is in a good agreement with the disk truncation radius 
inferred by other measurements, e.g. derived from the variability of the reflected 
component ([H], Sect. [7]). 

The right panel in Fig. [S| also demonstrates that there is considerable power at 
very high frequencies, at the level of uP^, ~ 1 — 3% at / ~ 10^ Hz. The millisecond 
range includes a number of important characteristic time scales in the vicinity of 
the compact object. The Keplerian frequency at the last marginally stable orbit for 
a non-rotating black hole is /k ~ 220 (AZ/IOM©)"^ Hz. The relativistic precession 
frequencies associated with Keplerian motion in the vicinity of the compact object 
are also in this range [54]. The typical value of the sound crossing time in the gas 
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pressure supported thin disk corresponds to the frequency of a few hundred Hz. The 
light crossing time for the region of ~ lOr^ corresponds to the frequency of ~ 1 
kHz. However, the observed power spectra break somewhere between / ~ 20 Hz 
(soft state) and / ~ 50 — 100 Hz (hard state) and do not exhibit any detectable 
features beyond these frequencies. The upper limit on the fractional rms of aperiodic 
variability in the 500-1000 Hz range is ~ 1%. Similar upper limit can be placed on 
the rms amplitude of narrow features anywhere in this range, ~ 0.9% assuming width 
V I Av = 20 (all upper limits at 95% confidence) [4^. This suggests that none of the 
above processes leads to a significant variability and, in particular, no resonances 
are present at their characteristic time scales. The nature of the high frequency 
breaks is unclear. It may be plausible to associate them with the properties of the 
emission mechanism of the hard component, for example with the details and time 
scales of heating and cooling of electrons in the Comptonization region. Obviously 
these properties are different for the hard component in the soft and hard state as 
indicated by the difference in the high frequency power density spectra shown in the 
right panel in Fig. [H] 

6.4 (Lack of) Variability in the disk emission 

The energy spectra of black holes in the soft state often have two spectral compo- 
nents, a soft blackbody-like component emitted by the optically thick accretion disk 
and a hard power law component of thermal or non-thermal origin formed in the 
optically thin media (Sect. |31 Figs. Ill3p . The study of the variability properties of 
these two components in black hole X-ray novae in the early 90's with the Ginga and 
GRANAT observatories led to the conclusion that most of the variability of the X- 
ray emission in the soft state is usually associated with the hard spectral component 
(e.g. |37II59|[T6| ). the thermal disk emission generally being much more stable. This 
conclusion has been later confirmed by a more rigorous and quantitative analysis 
[S] 1191 121j whose main results are summarized below. 

The energy dependence of the variability can be characterized by the frequency- 
resolved energy spectrum [45] • It is defined as a set of Fourier- amplitudes computed 
from light curves in different energy channels. The Fourier-amplitudes are integrated 
over the frequency range of interest and expressed in the units of flux. Its advantages 
over simple energy-dependent fractional rms is the possibility to use conventional 
(i.e. response folded) spectral approximations and to compare its shape with shape 
of various spectral components present in the average energy spectrum of the source. 
However it can not be always regarded as the energy spectrum of the variable com- 
ponent. For this to be possible, certain conditions must be fulflUed [21) . for example 
(i) independence of the shape of the frequency resolved spectrum on the Fourier- 
frequency and (ii) absence of time lags between flux variations at different energies. 
If these two conditions are satisfied, the spectral variability can be represented as 
F{E,t) = So{E)+I{t)S{E). 

These conditions are fulfilled in the soft state of Cyg X-1 (Fig. I^}: therefore the 
Fourier- frequency resolved spectrum represents the energy spectrum of the variable 
part of the X-ray emission. Remarkably, it coincides with the average source spec- 
trum at _E > 7 — 8 keV where the contribution of the disk emission becomes small 
fFig. nop. Moreover, source light curves in different energy channels allow a linear 
decomposition in the form F{E,t) = A{E) + I{t)B{E) 8 . The constant part of the 
source emission A{E) coincides with the spectrum of the accretion disk, while the 
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Fig. 9. Frequency resolved spectra of Cyg X-1 in the soft state (June 16-18, 1996). 
Tlie numbers in tlie upper-riglit corner of eacii panel indicate the median frequency. 
In each panel, the upper histogram shows the average spectrum, the lower histogram 
shows the frequency resolved spectrum at low frequencies, 0.002-0.033 Hz (the same 
data as in the left panel of Fig. I10|l . The spectra are plotted as ratio to a power law 
spectrum with photon index F — 2.5 and low energy absorption Nh ~ 6- 10^^ cm~^. 



spectrum of the variable part B[E) coincides with the frequency- resolved spectrum 
and can be described by the Comptonization model (Fig. llOp . Thus, on time scales 
from ~ 100 msec to ~ 500 sec (at least) the source variability is due to variations 
in the flux of the hard component, whose shape is kept constant in the course of 
these variations, the amplitude of variability of the disk emission being significantly 
smaller. 

This can be explained in the model of inward propagating fluctuations outfined 
in Sect. 16.11 The hard component is associated with the optically thin hot flow 
whose aspect ratio must be comparable to unity, h/r ~ 0.1 — 1. Therefore the vis- 
cous time scale in the hot flow is ~ 10^ — 10* times shorter than in the geometrically 
thin disk from which the soft component originates (cf. Sect. W^ . This makes the 
geometrically thick flow more "transparent" for the high frequency perturbations. 
Depending on the particular value of h/r and of the a-parameter, the viscosity or M 
fluctuations on the thermal and even dynamical time scales will be propagated in- 
ward without significant damping and will modulate the accretion rate at aU smaller 
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Fig. 10. Left: Average and frequency-resolved spectrum of Cyg X-1 in the 1996 
soft state. Right: Spectra of "constant" (open circles - A{E)) and "variable" (solid 
circles - B{E)) components, derived from the linear fits of the correlation between 
count rate in different channels. The normalization of the "variable" component 
B{E) is arbitrary. For comparison, the light grey (yellow in the color version) curve 
shows the spectrum of a multicolor disk emission with a characteristic temperature 
of 0.5 keV. The two upper spectra (solid squares) were averaged over the periods 
of time when the count rate above 9 keV was high and low respectively. The dark 
grey (green) lines show that these spectra can be reasonably well (within 10-15%) 
approximated by a model M{E) — A{E) + I * B{E) consisting of the stable and 
variable spectral components where / (the normalization of the variable component) 
is the only free parameter. The right panel is adopted from [8]. 



radii, including the region where X-ray emission is formed, leading to a significant 
modulation of X-ray flux. This is not the case for the geometrically thin disk where 
the viscous time scale is ~ 10'^ — lO'' times longer than the dynamical and thermal 
timescales; therefore high frequency perturbations will be damped (Sect. I6.1|l . It 
is plausible to expect that perturbations on the thermal and dynamical time scales 
have larger amplitude than perturbations on the viscous time scales, thus explaining 
the significantly smaller variability amplitude of the disk emission. 

The lack of variability of the disk emission and its interpretation are consistent 
with the conclusion of the Sect. 16.21 made from completely independent arguments. 
Namely, based on the location of the low-frequency break in the power density 
spectra of neutron stars we concluded that bulk of variability seen at low frequencies 
originates in the optically thin coronal flow with an aspect ratio of h/r ~ 0.1, rather 
than in the underlying geometrically thin disk. 

Finally, it should be noted that variations of the soft component may also arise 
from the variations of the disk truncation radius. Such variations may be absent 
if the geometrically thin disk extends to the last marginally stable orbit (3rg for 
a Schwarzschild black hole) but appear when the disk truncation radius is larger. 
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Indeed, detailed study of the 1996 soft state of Cyg X-1 used as an example here 
demonstrated that on a number of occasions, mostly at the beginning and in the end 
of soft state episodes, the soft component was strongly variable. Notably, the power 
density spectrum during these periods had a complex shape, significantly different 
form the simple power law shown in Figs l6l8l 

7 Variability of the reflected emission 

The reflected component (Sect.[4| arises from reprocessing of the Comptonized emis- 
sion in the accretion disk (Fig. therefore it should be expected to show some 
degree of variability, as the Comptonized radiation is strongly variable. The charac- 
teristic times of absorption/emission processes in the accretion disk are negligibly 
small, therefore the main factors defining the response of the reflected fiux to vari- 
ations of the Comptonized radiation are related to the light travel times (Fig. . 
Namely, they are: (i) the finite light travel time from the source of primary radiation 
to the reflector, i.e. from the Comptonization region to the accretion disk and (ii) 
the finite size of the reflector. The flrst will introduce a time delay between vari- 
ations of the reflected and Comptonized components. The amplitude of this delay 
is Td ~ Td/c ~ 10 (ra/WOrg) (M/10Mq) msec. The finite size of the disk itself will 
lead to suppression of the high frequency variations in the reflected component - 
the accretion disk acts as a low-pass filter. It seems to be possible to estimate the 
cut-off frequency from the size of the accretion disk making the main contribution 
to the refiected fiux, Avd ~ rd, thus leading to f^ut ~ ^^d^ ~ {Ard/c)~^ ~ 100 Hz, 
i.e. beyond the frequency range of the bulk of observed variability (Fig.[8)|. However, 
calculations of the transfer function [T5] show that at the frequency fcut ~ '^^d^ 
the variability signal is suppressed by a significant factor, of ~ 10 — 20, whereas a 
noticeable suppression of variability, by a factor of ~ 2 or more occurs at frequen- 
cies ~ 10 times lower. As these effects directly depend on the mutual location and 
geometry of the Comptonization region and accretion disk, their observation is a 
powerful tool in studying the geometry of the accretion flow. 

As suggested in [451 119) , the variability of the refiected emission can be studied 
using the methods of frequency-resolved spectroscopy (Sect. 16. 4p . This method is 
based on the analysis of Fourier amplitudes and therefore washes out phase infor- 
mation, precluding the study of time delay effects. However, it provides a convenient 
way to explore the frequency dependence of the variability of the iron line flux and 
to compare it with that of the continuum emission which is dominated by the Comp- 
tonized radiation. Using this method, [451 119) showed that in the soft spectral state, 
variations of the reflected component have the same frequency dependence of the 
rms amplitude as the Comptonized emission up to frequencies ~ 30 Hz (Fig. I12|) . 
This would be expected if, for instance, the reflected flux was reproducing, with a 
flat response, the variations of the primary radiation down to time-scales of ~ 30 — 50 
ms. The sensitivity of their analysis was insufficient to study shorter time-scales. In 
the hard spectral state, on the contrary, the variability of the reflected fiux is sig- 
nificantly suppressed in comparison with the direct emission on time-scales shorter 
than ~ 0.5 — Is. These findings are to be compared with the predictions of the 
simple model for the time response of the disk to variations of the primary emission. 
Assuming that suppression of the short-term variability of the reflected emission is 
caused by the flnite light-crossing time of the accretion disk, one can estimate the 
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Fig. 11. The two effects defining the response of the reflected emission to variations 
of the Comptonized emission are the flnite light travel time from the hot inner 
flow to the accretion disk and the finite size of the accretion disk itself, Ar^/c ~ At^ 
(left panel). While the first factor introduces a time delay between variations of the 
reflected and Comptonized components, the finite size of the disk itself leads to the 
suppression of high frequency variations in the reflected emission. The right panel 
shows the Green function of the time response of the geometrically thin disk around 
a IOMq black hole (from [H]). The numbers at the curves mark the disk truncation 
radius and inclination. 

truncation radius of the accretion disc, Vd ~ lOOrg in the hard spectral state and 
rd < lOr^ in the soft spectral state (Fig:.ll2p. This agrees well with the interpretation 
of the spectral states in black holes in the "sombrero" geometry of the accretion flow 
(Sect. El). 

In order to study the time delay in the reflected emission, one would need to 
separate the reflected component from the main emission in time-resolved energy 
spectra. This can be done most easily for the fluorescent iron line, while the reflection 
continuum is more difficult to separate due to its breadth in the energy domain. An 
attempt to perform such an analysis based on RXTE/PCA data of Cyg X-1 in 
the hard state is presented in the right panel in Fig. 1121 In each bin of the light 
curve with 16 msec resolution, the spectrum in the 3-20 keV band was linearly 
decomposed in to power law component, refiected continuum and 6.4 keV iron line 
emission. Shown in Fig. [12] is the cross-correlation of the light curves of the iron 
line flux and the power law component. The cross-correlation function shows a clear 
asymmetry, suggesting a time delay of the order of ~ 10— 15 msec. The amplitude of 
the possible time delay corresponds to a disk truncation radius of ~ lOOr^, in good 
agreement with the number obtained from the analysis of high-frequency variations 
in the iron line flux (left panel in Fig. I12|) and also with the numbers tentatively 
suggested by spectral analysis. 

The results shown in Fig. [T^] seem to suggest a rather consistent picture and to 
favor to the "sombrero" -type configuration of the accretion flow, with the spectral 
state transition being related to a change of the disk truncation radius. However, 
a caveat is in order. While the rms amplitude behavior of the iron line emission 
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Fig. 12. Suppression of high frequency variations (left) and time delay (right) in the 
reflected emission. The left panel shows the equivalent width of the iron fluorescent 
line {— ratio of the fractional rms of the iron line flux to that of the underlying 
continuum) vs. Fourier frequency. The model curves are for an isotropic point source 
at the height h = lOvg on the a:xis of a flat disk with inner truncation radius of 10, 100 
and lOOOrg (assuming a IOM0 black hole) and with inclination angle of 50 degrees. A 
narrow line was assumed in calculations. The right panel shows the cross-correlation 
function of the iron line emission and the power law continuum. Positive lags mean 
delay of the line emission with respect to the continuum. 

is a rather robust observational result, the search of the time delay of the reflected 
component requires separation of line and continuum emission components. On ~ 10 
msec time scales, this can not be done through direct spectral fltting because of 
insufficient statistics, even with the large collecting area of the PCA instrument 
aboard RXTE, and requires more sophisticated data analysis techniques, for example 
the one used to produce the cross-correlation function shown in Fig. [12] Secondly, 
although the interpretation of these results in terms of finite light travel times is the 
most simple and straightforward, alternative scenarios are also possible as discussed 
in length in |19) . Nevertheless, the former interpretation is, in my view, the simplest 
and most attractive one and is further supported by the results of the spectral 
analysis, as described in the next section. 



8 R — r and other correlations 

Observations show that spectral and timing parameters of accreting black holes often 
change in a correlated way, e.g. |17ll43lHD] . One of the most significant correlations 
is the one between the photon index of the Comptonized spectrum, the amplitude 
of the reflected component and the characteristic frequencies of aperiodic variability 
fFigs. [T3ll6p . The correlation between spectral slope and reflection amplitude is also 
known as R — F correlation |17l I67j . Its importance is further amplified by the fact 
that it is also valid for supermassive black holes (Fig. [T4)l [68) . 
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Fig. 13. The R — F correlations between the photon index F of the Comptonized 
radiation and the relative amplitude of the reflected component R in three "well- 
behaved" black hole systems Cyg X-1, GX 339-4 and GS 1354-644. The solid and 
dashed lines show the the dependence r{R) expected in the disk-spheroid and in 
the plasma ejection models discussed in the text. Adopted from |18] 

8.1 R — r correlation 

The spectrum formed by the unsaturated Comptonization of low frequency seed 
photons with characteristic temperature Ttb on hot electrons with temperature Te 
has a nearly power law shape in the energy range from ~ 3kTi,t to ~ kT^ [56j . 
For the parameters typical for black hole X-ray binaries in the hard spectral state, 
this corresponds to the energy range from ~ 0.5 — 1 keV to ~ 50 — 100 keV (e.g. 
Fig. [21). The photon index F of the Comptonized spectrum depends in a rather 
complicated way on the parameters of the Comptonizing media, primarily on the 
electron temperature and the Thompson optical depth [56]. It is more meaningful to 
relate F to the Comptonization parameter y or, nearly equivalently, to the Compton 
amplification factor A. The latter describes the energy balance in the corona and is 
defined as the ratio of the energy deposition rate into hot electrons and the energy 
flux brought into the Comptonization region by soft seed photons. The concrete 
shape of the F{A) relation depends on the ratio Tbb/T^ of the temperatures of the 
seed photons and the electrons, the Thomson optical depth and the geometry, but 
broadly speaking, the higher the Compton ampliflcation factor, the harder is the 
Comptonized spectrum 57 , 1111 1231 116j . This is illustrated by the results of Monte- 
Carlo simulations shown in Fig. [TJl 

The strength of the reflected component in the spectrum depends on the fraction 
of the Comptonized radiation intercepted by the accretion disk (Sect. [4)1. The latter 
is deflned by the geometry of the accretion flow, namely, by the solid angle fidisk 
subtended by the accretion disk as seen from the corona. In addition, the spectrum 
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of the reflected emission depends on the ionization state of the disk, in particular its 
low energy part which is formed by the interplay between Thomson scattering and 
photoabsorption and fluorescence by metals. The problem is further complicated by 
the fact that the ionization state of the disk can be modified by the Comptonized 
radiation. 

Observations show that there is a clear correlation between the photon index of 
the Comptonized radiation F (i.e. the Comptonization parameter) and the relative 
amplitude of the reflected component R (Fig. I13|) . Softer spectra (lower value of 
the Comptonization parameter y and of the Compton amplification factor A) have 
stronger refiected component, revealing itself, for example, via a larger equivalent 
width of the iron fiuorescent line. The existence of this correlation suggests that 
there is a positive correlation between the fraction of the Comptonized radiation 
intercepted by the accretion disk and the energy fiux of the soft seed photons to the 
Comptonization region [67]. This is a strong argument in favor of the accretion disk 
being the primary source of soft seed photons to the Comptonization region. Indeed, 
in the absence of strong beaming effects a correlation between Odisk and the seed 
photons flux should be expected since an increase of the solid angle of the disk seen 
by the hot electrons {— Qdisk) should generally lead to the increase of the fraction 
of the disk emission reaching the Comptonization region. 



8.2 Toy models 

We illustrate the above considerations with two simple and idealized models having 
a different cause of change of fldisk- In the flrst, the disk-spheroid model (cf. "som- 
brero" conflguration. Sect. [2]), an optically-thin uniform hot sphere with radius Vsph, 
the source of the hard Comptonized radiation is surrounded by an optically-thick 
cold disk with an inner radius rdisk, ^disk depending on the ratio rdisk/rgph- Propa- 
gation of the the disk towards/inwards the hot sphere (decrease of Tdisk/raph) leads 
to an increase of the reflection scaling factor R, a decrease of the the Compton am- 
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Fig. 15. Relation between tlie plioton index of Comptonized radiation 7^ and the 
Compton amplification factor A. The symbols show results of Monte-Carlo simula- 
tions assuming spherical geometry for different values of parameters of the Comp- 
tonization region and soft seed photons. The solid lines are calculated using eq.([l}. 
Adopted from ^ISj 



plification factor A and a steepening of the Comptonized spectrum. In such context 
the model was first studied by Zdziarski, Lubinski & Smith [67] . In the second, the 
plasma ejection model, the value of the fidiak is defined by the intrinsic properties 
of the emitting hot plasma, particularly by its bulk motion with mildly relativistic 
velocity towards or away from the disk, which itself remains unchanged [6j. In the 
case of an infinite disk, values of the refiection scaling factor R below and above 
unity correspond to the hot plasma moving respectively away from and towards the 
disk. 

Both models predict a relation between refiection R and Compton amplification 
factor A which can be translated to r{R) given a dependence r{A) of the photon 
index of the Comptonized spectrum on the amplification factor. The relation between 
7^ and A can be approximated by: 

This formula is based on a representation of the Comptonized spectrum by a power 
law in the energy range kTbb ~ kTe and takes into account that a fraction e~^^ 
of the soft radiation will leave the Comptonization region unmodified. Despite its 
simplicity, it agrees with the results of the Monte-Carlo calculations with reasonable 
accuracy for optical depth tt ~ 1 and Ttt/Te ~ 10~^ — 10^"* (Fig. [T^ . 

The expected r{R) relations are shown in Fig. 1131 With a proper tuning of the 
parameters, both models can reproduce the observed shape of the r{R) dependence 
and in this respect are virtually indistinguishable. The models plotted in Fig. [T3] 
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were calculated with the following parameters: the disk-spheroid model assumes 
the disk albedo a — 0.1, Thomson optical depth of the cloud tt ~ 1 and the ratio 
of the temperature of the seed photons to the electron temperature Tbb/Te = 10~*. 
I note that the latter value is too small for stellar mass black holes, a more realistic 
one being in the range Tbb/Te — 10~^. However, this does not invalidate the model 
as a number of significant effect are ignored in this calculation which can modify the 
r{R) dependence, for example disk inclination, gravitational energy release in the 
disk etc. The exact importance of these factors is yet to be determined. The plasma 
ejection model parameters are: a = 0.15, tt = 1, Tbb/Te — 3 • 10~^ and fis ~ 0.3. 
The observed range of the reflection R ~ 0.3 — 1 and the slope P ~ 1.5 — 2.2 can 
be explained assuming variation of the disk radius from rdisk ~ '"aph to rdisk ~ in 
the disk-spheroid model or variation of the bulk motion velocity from v ~ 0.4c away 
from the disk to ~ in the plasma ejection model. 

Of course these models are very simple and schematic and the real configuration 
of the accretion flow is likely to be far more complex. They are presented here with 
the only purpose to demonstrate that simple geometrical considerations can suc- 
cessfully explain the observed correlation between parameters of the Comptonized 
and reflected emission in black holes. More sophisticated scenarios are considered 
for example in [29] . 

8.3 Characteristic frequencies of variability 

As discussed in the Sect. 16.31 the power density spectra of accreting black holes above 
/ > lO"'^ Hz have a number of bumps and peaks which deflne several frequencies 
characterizing the variability time scales in the accretion flow. These frequencies usu- 
ally correlate with each other therefore almost any of them may be used to represent 
characteristic variability timescales. As illustrated by Fig. 1161 a tight correlation ex- 
ists between the reflection amplitude R and the characteristic variability frequencies 
- an increase of the amplitude of the reflected component in the energy spectrum is 
accompanied by an increase of the variability frequencies. This correlation covers a 
remarkably broad dynamical range, nearly two orders of magnitude in frequency. 

Although the precise nature of the characteristic noise frequencies is still un- 
known, it is plausible that they are associated with the Keplerian and viscous time 
scales of the disk and corona at various characteristic radii, for example, at the 
truncation radius of the disk (Sect. [6^ . If this is true, the correlation between R 
and V can be easily understood, at least qualitatively. Indeed, in the truncated disk 
picture the increase of reflection is caused by the inward propagation of the inner 
disk boundary, hence, it is accompanied by an increase of the Keplerian frequency 
at the disk truncation radius and corresponding increase of the characteristic noise 
frequencies. 

8.4 Doppler broadening of the iron line 

The spectrum of the emission, reflected from a Keplerian accretion disk, is modified 
by special and general relativity effects [T3] , in particular the width of the fiuorescent 
line of iron is affected by the Doppler effect due to Keplerian motion in the disk (Sect. 
[4} . If the increase of the reflection amplitude is caused by the decrease of the inner 
radius of the accretion disk, a correlation should be expected between the amplitude 
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Fig. 16. The correlations between the amplitude of the reflected component R and 
characteristic frequencies of aperiodic variability. The frequency of the second peak 
in the u Pi, plot (Fig[8} was used to represent the latter. 



of reflected emission and its Doppler broadening, in particular the Doppler width of 
the iron line. Such a correlation is a generic prediction of the truncated disk models 
and might be used to discriminate between different geometries of the accretion flow. 
The energy resolution of RXTE/PCA, whose data has been used for this study [20] . 
is not entirely adequate for the task to accurately measure the relativistic smearing 
of the reflection features. However, the data shown in Fig. 1171 suggests a correlated 
behavior of the reflection and the Doppler broadening of the fluorescent line of iron. 

Speculating further, if the characteristic frequencies of variability are propor- 
tional to the Keplerian frequency at the inner boundary of the disk, they should 
scale as 

UQPo cxujK cc r~fj^ 

As the reflected emission is likely to originate primarily from the innermost parts of 
the accretion disk, closest to the source of Comptonized radiation, the effect of the 
Doppler broadening should be proportional to the Keplerian linear velocity at the 
inner edge of the disk: 

Vk . . -1/2 . . 

—pr oc Sim oc r^. ',, smi 

E c 

Therefore, one might expect that the characteristic frequencies of variability and the 
Doppler broadening of the fluorescent line should be related via: 

(AE/E? 

VQPO OC — r-T- 

Mbh sm^'z 



where Mbh is the black hole mass and i is inclination of the binary system. The 
PCA data indicate that such dependence might indeed be the case (Fig. 1171 right 
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Fig. 17. Relation between the Doppler broadening of the iron hne and the reflec- 
tion scahng factor (left) and the characteristic frequency of variability (right). The 
straight lines in the right panel show the dependence AE/E oc i^Qpo expected in 
the truncated disk picture if the characteristic frequencies of variability were pro- 
portional to the Keplerian frequency at the inner boundary of the accretion disk. 

panel). However, an independent confirmation by observations with higher energy 
resolution instruments is still needed. 

Thus, observations speak in favor of the truncated disk scenario. However this 
point of view is not universally accepted and alternative interpretations and scenar- 
ios are being investigated. The counter arguments are based, for example, on the 
possible detection of the relativistic broad Fe K line and cool disk emission compo- 
nent in the hard state spectra of several black holes, suggesting that the optically 
thick disk might be present in the vicinity of the compact object in the hard state 
as well. For the detailed discussion the interested reader is referred to the original 
work, e.g. ^5] and references therein. 



9 Comparison with neutron star binaries 

Neutron star radii are most likely in the range ~ 10 — 15 km, i.e. of the order of ~ 
3 — irg for a ~ IAMq object. This is comparable to the radius of the last marginally 
stable Keplerian orbit around a non-rotating black hole, hence the efficiency r] = 
Lx jM(? of accretion onto a neutron star is not much different from that on a black 
hole. Therefore at comparable M accreting black holes and neutron stars would 
have comparable X-ray luminosities (but see below regarding the contribution of the 
boundary layer). The Eddington luminosity limit, however, is proportional to the 
mass of the central objects, and the typical maximum luminosity is by Mbh /Mns ~ 
5—10 times smaller for neutron stars. This is in general agreement with observations 
of peak luminosities of black hole transients and their comparison to luminosities of 
persistent and transient neutron star systems in the Milky Way. 
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Fig. 18. The left panels show the geometry of the spreading layer on the surface of 
the neutron star (upper) and radial dependence of the angular velocity of accreting 
matter (lower). Adopted from [25]. On the right: Hard state spectra of several weakly 
magnetized accreting neutron stars and of a black hole Cyg X-1. Based on the data 
of RXTE observations. 

The fact that the size of a neutron star is of the order of ~ Srg, i.e. is comparable 
to the radius of the last marginally stable Keplerian orbit around a black hole, also 
suggests that the structure of the accretion disk in both cases may be similar. This 
is indirectly confirmed by the existence of two spectral states in accreting neutron 
stars, which properties are qualitatively similar to black holes (Fig. [!}. A more 
direct argument is presented by the similarity of the accretion disk spectra in the 
soft state, which can be described by the same spectral models in black holes and 
neutron stars (e.g.^2lj'). Moreover, the power density spectra of the accretion disk 
in both types of systems follow a Pi^ oc power law with a close to ~ 1 (Sect. 
Eland [21]). The strong magnetic field which may exist around young neutron stars 
can change the picture, causing disruption of the accretion disk at a large distance 
from the neutron star and modifying dramatically the structure of the accretion flow 
inside the magnetospheric radius Vm 3> rjvs- This may lead to the phenomenon of 
X-ray pulsations common in high-mass X-ray binaries and is not considered in this 
chapter. 

The main qualitative difference between the two types of compact objects is 
obviously the existence of a solid surface of the neutron star which is absent in 
the case of a black hole. Neutron star rotation frequencies are typically in the few 
hundred Hz range, i.e. ~ several times smaller than the Keplerian frequency near 
its surface, vk ~ kHz. Therefore a boundary or spreading layer will appear near 
the surface of the neutron star where accreting matter decelerates from Keplerian 
rotation in the accretion disk down to the neutron star spin frequency and settles 
onto its surface (Fig.[T8l left panels). In Newtonian approximation, half of the energy 
of the test particle on a Keplerian orbit is in the form of kinetic energy of the 
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Keplerian rotation, hence for a non-rotating neutron star half of the energy release 
due to accretion would take place in the boundary /spreading layer. The effects of 
general relativity can increase this fraction, e.g. up to ~ 2/3 in the case of a neutron 
star with radius rjvs = 3rg [531 151) . Rotation of the neutron star and deviations of 
the space-time geometry from Schwarzschild metric further modify the fraction of 
the energy released on the star's surface. A luminous spectral component emitted 
by the boundary layer will exist in the X-ray spectrum of an accreting neutron star, 
in addition to the emission from the accretion disk. The shape of the spectra of 
luminous neutron stars (soft state spectrum in Fig. [1} suggests that both accretion 
disk and boundary layer are in the optically thick regime. As the luminosities of 
the two components are comparable, but the emitting area of the boundary layer is 
smaller than that of the accretion disk, its spectrum is expected to be harder. 

Looking from a different angle we may say that because of the presence of the 
stellar surface, at the same M neutron stars are approximately twice more luminous 
than black holes. Indeed, in the case of the black hole accretion the kinetic energy of 
Keplerian motion at the inner edge of the accretion disk is nearly all advected into 
the black hole. If the central object is a neutron star, this energy is released in the 
boundary /spreading layer on its surface, approximately doubling the luminosity of 
the source|j An interesting consequence of this is that the ratio of LEdd/M is larger 
for neutron stars than for black holes. 

Due to the similarity of the spectra of the accretion disk and boundary layer, 
the total spectrum has a smooth curved shape, which is difficult to decompose into 
separate spectral components. This complicates analysis and interpretation of the 
neutron star spectra and, in spite of a very significant increase in the sensitivity of 
X-ray instruments made in the last decade, still often leads to ambiguous and contra- 
dicting results, even based on physically motivated spectral models. This degeneracy 
may be removed with the help of timing information. It has been noticed already in 
the early 80's that the variability patterns may be different for the boundary layer 
and accretion disk ^36,. Further progress has been made almost twenty years later, 
thanks to the large collecting area of the PGA instrument aboard RXTE and the use 
of novel data analysis techniques. Using the method of frequency-resolved spectral 
analysis it has been shown that aperiodic and quasi-periodic variability of bright 
LMXBs - atoll and Z- sources, on ~ sec - msec time scales is caused primarily by 
variations of the luminosity of the boundary layer [5T]. It was also shown that the 
boundary layer spectrum remains nearly constant in the course of the luminosity 
variations and is represented to certain accuracy by the Fourier frequency resolved 
spectrum. This permits to resolve the degeneracy in the spectral fitting and to sep- 
arate contributions of the boundary layer and disk emission (Fig. 1191 left panel). 
Interestingly, the spectrum of the boundary layer emission has the same shape in 
different objects and is nearly independent upon the global mass accretion rate in 
the investigated range of M ~ (0.1 — l)A4"Edd and in the limit of M ~ AfEdd is 
close to a Wien spectrum with kT ~ 2.4 keV (Fig. 1191 right panel). Its indepen- 
dence on the global value of M lends support to the theoretical suggestion by [25] 



^ A more accurate consideration should take into account geometry of the problem, 
namely the emission diagrams and orientation with respect to the observer of 
the emitting surfaces of the boundary layer and accretion disk. This makes the 
boundary layer contribution to the observed emission dependent on the inclination 
of the binary system, see for example [21] . 
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Fig. 19. Left: Points with error bars showing total and frequency resolved-spectra 
of GX 340+0. The dashed (blue in the color version of the plot) and upper solid 
(red) histograms show the disk and boundary layer spectra, the latter computed as 
a difference between the (observed) total and (predicted) accretion disk spectrum. 
The lower solid histogram is obtained from the upper scaling it to the total energy 
flux of the frequency-resolved spectrum. Adopted from 21 . Right: Fourier-frequency 
resolved spectra, corrected for the interstellar absorption (^boundary layer spectra) 
of 5 luminous accreting neutron stars. For 4U 1608-52, the frequency-resolved spec- 
trum of the lower kHz QPO is shown. All spectra were corrected for the interstellar 
absorption. The thick short-dashed line shows the best fit Comptonization model 
with kTs = 1.5, fcTe = 3.3 keV, r = 5. The thin long-dashed line shows a blackbody 
spectrum with temperature feTbb = 2.4 keV. Adopted from [46] . 

that the boundary layer is radiation pressure supported. With this assumption, one 
can attempt to measure gravity on the neutron star surface and hence M/R^ , from 
the shape of the boundary layer spectrum, similarly to the photospheric expansion 
(i.e. Eddington limited) X-ray bursts. This gives results within the range of values 
obtained by other methods [461 152j . 

As neutron stars have smaller mass than black holes, the linear scale corre- 
sponding to the gravitational radius is smaller, Vg oc M. This has two important 
consequences. Firstly, the surface area of the emitting region is oc M^, therefore 
at the same luminosity the temperature of black body emission will be larger, 
Tbb oc M~^^'^ . The spectrum is further modified by the contribution of the emis- 
sion from the boundary layer, which has yet smaller emitting region and harder 
spectrum. This is illustrated by Fig. [1]- the soft state spectrum of the neutron star 
4U 1705-44 is noticeably harder than that of the black hole Cyg X-1. 

Secondly, the smaller linear scale in neutron star accretion shifts the character- 
istic frequency scales by a factor oc M~^, suggesting ~ 5 times higher frequencies 
of variability in neutron star systems. This is confirmed by observations as demon- 
strated in the extensive comparison of black hole and neutron star power density 
spectra in [60]. In addition, characteristic time scales in the boundary or spreading 
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layer are significantly shorter than those in the accretion disk and are in the tens 
of kHz frequency domain. It has been suggested that a very high frequency com- 
ponent may exist in the power density spectra of neutron stars associated with the 
turbulence in the spreading layer [BOlUS]. If detected, this may become a unique 
diagnostics tool of physical conditions in the spreading layer on the surface of the 
neutron star. So far, only upper limits have been obtained, they are at the level of 
~ 10^*^ fractional rms [60] . 

The relatively cold surface of the neutron star is a source of copious soft pho- 
tons. This is mostly relevant in the hard spectral state - the low energy photons 
emitted by the neutron star surface result in a more efficient cooling of electrons in 
the Comptonization region, thus changing its energy balance and Comptonization 
parameter (luminosity enhancement factor, cf. Sect. I8.l( l. Consequently, the energy 
spectra of neutron stars in the hard state are significantly softer than those of black 
holes - they have large spectral index (smaller Comptonization parameter) and a 
smaller value of the high energy cut-off (lower electron temperature). Both these 
effects are clearly seen in the spectra shown in Figs. I1I18I 

The existence of soft and hard spectral states in neutron stars suggests that, 
similarly to black holes, the accretion flow can change its configuration from optically 
thin to optically thick. A remarkable fact is that this change seems to occur to both 
accretion disk and boundary layer (quasi) simultaneously. Indeed, to my knowledge, 
no two-component (soft -I- hard) spectrum has been observed in the case of accreting 
neutron stars so far. 
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